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In this talk, we introduce our recently completed next-to-leading order (NLO) global analysis of 
the nuclear parton distribution functions (nPDFs) called EPS09 — a higher order successor to 
the well-known leading-order (LO) analysis EKS98 and also to our previous LO work EPS08. 
As an extension to similar global analyses carried out by other groups, we complement the data 
from deep inelastic / +A scattering and Drell-Yan dilepton measurements in pH-A collisions by 
inclusive midrapidity pion production data from dn-Au collisions at RHIC, which results in better 
constrained gluon distributions than before. The most important new ingredient, however, is the 
detailed error analysis, which employs the Hessian method and which allows us to map out the 
parameter- space vicinity of the best-fit to a collection of nPDF error sets. These error sets provide 
the end-user a way to compute how the PDF-uncertainties will propagate into the cross sections 
of his/her interest. The EPS09 package to be released soon, will contain both the NLO and LO 
results for the best fits and the uncertainty sets. 



4th international workshop High-pT physics at LHC 09 
February 4-7, 2009 
Prague, Czech Republic 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



EPS09 - Global NLO analysis of nuclear PDFs and their uncertainties 



H. Paukkunen 



1. Introduction 

The global analyses of the free-nucleon parton distribution functions (PDFs) are based on three 
fundamental aspects of QCD: asymptotic freedom, coUinear factorization and scale evolution of the 
PDFs. These properties allow to compute the inclusive hard process cross-sections, schematically 

OAB^h+X = (G') ® Oij^h+X ® ff{Q^) + ^ (e') " , 

'J 

where fis denote the universal, process-independent PDFs obeying the DGLAP [|l|] evolution equa- 
tions, and Oij^h+x are perturbatively computable coefficients. This approach has proven to work 
extremely well with increasingly more different types of data included in the analyses. In the 
case of bound nucleons the validity of factorization is not as well established (see e.g. Ref. 
but it has nevertheless turned out to provide a very good description of the world data from deep 
inelastic scattering (DIS) and Drell-Yan (DY) dilepton measurements involving nuclear targets 
[^, ^ ^ 1^ 0, ^ : in the nuclear environment the shape of the PDFs, however, is different from the 
free-nucleon PDFs. Here, we give an overview of our recent NLO global analysis of the nuclear 
PDFs (nPDFs) and their uncertainties [^. 

2. Analysis method and framework 

Our analysis follows the usual global DGLAP procedure: 

• The PDFs are parametrized at a chosen initial scale imposing the sum rules. In 

this work we do not parametrize the absolute nPDFs, but the nuclear modification factors 
Rf{x, Qq) encoding the relative difference to the free-proton PDFs through 

y;^(^,G^)^/?f(.,G^)/;"^Q^-^^(x,G2). (2.i) 

Above, ff^'^^^^-^^{x,Q^) refers to the CTEQ6.1M set of free proton PDFs [|l|] in the MS 
scheme and zero-mass variable flavour-number scheme. We consider three different modifi- 
cation factors: Ry{x,Ql) for both u and d valence quarks, R^{x,Qq) for all sea quarks, and 
Rq{x,QI) for gluons. 

• The absolute PDFs at the parametrization scale Qq are connected to other perturbative 
scales > Qq through the DGLAP evolution. In this way, also the nuclear modification 
factors Rf{x,Q^) become scale-dependent and the initial flavour-independence may also dis- 
appear. An efficient numerical solver for the parton evolution is an indispensable tool for 
any global analysis, but in the case of nPDFs this is even more critical as we always need 
to perform the evolution separately for 13 different nuclei — an order of magnitude more 
computation time than what is needed in a free-proton analysis. 

• The cross-sections are computed using the factorization theorem. In the present analysis 
the DIS and DY data constitute the main part of the experimental input, but we also em- 
ploy the midrapdity Ti^-production data measured in d-i-Au collisions at RHIC. Inclusion of 
the TT^-data provides important further constraints for the gluon modification that are partly 
complementary to the DIS and DY data. 
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The computed cross-sections are compared with the experimental data, and the initial 
parametrization is varied until the best agreement with the data is found. Our definition 
for the "best agreement" is based on minimizing the following ;f ^-function 

X\{a}) ^ IwA,z^({a}) (2.2) 



N 

2 'fNDi-Ti{{a}) 



/ ieN 



(2.3) 



Within each data set N, Dj denotes the experimental data value with a, point-to-point un- 
certainty, and Ti is the theory prediction corresponding to a parameter set {a}. If an over- 
all normalization uncertainty a™™ is specified by the experiment, the normalization factor 
/a? G [1 — CJ™™, 1 + CJ^™] is introduced. Its value is determined by minimizing Xn ^^'^ the 
final /n is thus an output of the analysis. The weight factors wn are used to amplify the 
importance of those data sets whose content is physically relevant, but contribution to X'^ 
would otherwize be too small to have an effect on the automated x^ minimization. 

In addition to finding the parameter set {a°} that optimally fits the experimental data, quantify- 
ing the uncertainties stemming from the experimental errors has become an increasingly important 
topic in the context of global PDF analyses. The Hessian method [ jlTl ] provides a practical way of 
treating this issue. It is based on a quadratic approximation for the x^ around its minimum Xq, 

X' ^ ^0 +1 ^^f^(«< - - = + LM^i - -a)). (2.4) 

which defines the Hessian matrix H. Non-zero off-diagonal elements in the Hessian matrix are a 
sign of correlations between the original fit parameters, invaUdating the standard (diagonal) error 
propagation formula 

(AX)2 ={^-8a^\(^- 5a2) \- (2.5) 
\oa\ J \oa2 J 

for a PDF-dependent quantity (cross section) X. Therefore, it is useful to diagonalize the Hessian 
matrix, such that 

X^^xl + Y^zj, (2.6) 

i 

where each z, is a certain linear combination of the original parameters around In these 

variables, the usual form of the error propagation 

dX ^ V fdX ^2 



stands on a much more solid ground. How to determine the size of the deviations 8zk is, however, 
a difficult issue where no universally agreed procedure exists. We consider 15 parameters in the 
expansion ( ^^ and the prescription which we adopt here is to choose each 5zk such that x^ grows 
by a certain fixed amount ^X^- Requiring each data set to remain close to its 90%-confidence 
range, we end up with a choice A;^^ = 50 (see Appendix A of Ref. [0] for details). 
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Much of the practicaUty of the Hessian method resides in constructing PDF error sets, which 
we denote by 5^. Each 5^ is obtained by displacing the fit parameters to the positive/negative 
direction along Zk such that grows by the chosen A;t^^ = 50. Approximating the derivatives in 



Eq. ( p.7| ) by a finite difference, we may then re-write the error formula as 

1 
4 



m^=\L[nst)-nsu)\\ (2.8) 
k 



where ^(5^) denotes the value of the quantity X computed by the set 5^. If the lower and upper 
uncertainties AX* differ, they should be computed separately, using the prescription [12] 

(AX+)2 « £ [max{Z(5+) -X{S''),X{S,) -X(5°),0}]' (2.9) 

k 

(AZ-)2 « £[max{X(5°)-X(5+),X(5°)-X(5,-),0}]', 

k 

where denotes the best fit. 

Along with the grids and interpolation routine for the best NLO and LO fits, our new release 



a computer routine called EPS09 [13] — contains also 30 eiTor sets required for computation of 



the uncertainties to nuclear cross-section ratios similar to those we present in the following section. 



However, as can be understood from the definition (gJJ), the total uncertainty in the absolute nPDFs 
is a combination of uncertainties from the baseline set CTEQ6.1M and those from the nuclear 
modifications Rf. Therefore, if an absolute cross-section is computed, also the CTEQ6.1M error- 
sets (with the EPS09 central set) should be included when calculating its uncertainty. 

3. Results 

In this section we present a selection of results from the NLO analysis. First, in Fig. |l| we 
plot the obtained modifications at two scales, at = 1.69GeV^ and at = lOOGeV^, which is 
to demonstrate the scale-dependence of the modifications and their uncertainties. One prominent 
feature to be noticed is that even if there is a rather large uncertainty band for the small-^; gluon 
modification R^ at Q^, the scale evolution tends to bring even a very strong gluon shadowing close 
to no shadowing at ~ lOOGeV^ — a very clear prediction of the DGLAP approach. 

The DIS data constitutes the bulk of the experimental data available for the global analysis of 
nPDFs. In Fig. ^ we show some of the measured nuclear modifications for (Z + A) DIS structure 
functions with respect to Deuterium, 

^^.(-,e^) = #|^, (3.1) 

and the comparison with our EPS09 NLO results. The shaded bands denote the uncertainty derived 
from the 30 EPS09NLO error sets and, as should be noticed, their width is comparable to the 
error bars in the experimental data. This a posteriori supports the validity of our procedure for 
determining the Ax^- Similar conclusion can be drawn upon inspecting the nuclear effects in the 
Drell-Yan data 

/?,,(.,„M)=^^^^,^^^^^^^^^, (3.2) 
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Figure 1: The nuclear modifications Rv, Rs, Rg for Lead at our initial scale Qq — 1.69GeV and at = 
lOOGeV^. The thick black lines indicate the best-fit results, whereas the dotted green curves denote the 
individual error sets. The shaded bands are computed using Eq. (2.9). 




Figure 2: The calculated NLO Rp^ (x, Q^) (filled symbols and error bands) compared with the NMC 95 [ p4t 
and the reanalysed NMC 95 [O] data (open symbols). 



where is the invariant mass of the lepton pair and xi 2 = \JM^ jse^^ {y is the pair rapidity). 
Comparison to the E772 and E866 data is shown in Fig. ||. Let us mention that the E772 data 
at X2 > 0. 1 carry some residual sensitivity also to the gluons and sea quarks, which has not been 
noticed before. 

The nuclear modification for inclusive pion production in d+Au collisions relative to p+p is 
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Figure 3: The computed NLO R^^{x,M^) (filled squares and error bands) as a function of xi and xi com- 
pared with the E866 [|l7|| and the E772 [|l6|] data (open squares). 
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Figure 4: The computed /?dAu (thick black line and blue error band) at y = for inclusive pion production 
compared with the PHENIX [|l8|] data (open squares). The error bars are the statistical uncertainties, and the 
yellow band indicates the point-to-point systematic errors. The PHENIX data have been multiplied by the 
optimized normalization factor = 1.03, which is an output of our analysis. The STAR data [ [l9| ] (open 
circles), multiplied by fy = 0.90, are also shown although they were not included in the EPS09 analysis. 



defined as 

(A^coii) d^Nl^/dpTdy d^af/dpTdy ' ^ ' ' 

where pT,y we. the transverse momentum and rapidity of the pion, and {N^ow) denotes the number 
of binary nucleon-nucleon collisions. A comparison with the PHENIX and STAR data is shown in 
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Fig. ^ and evidently, the shape of the spectrum gets very well reproduced by our parametrization. 
The downward trend toward large pT provides direct evidence for an EMC-effect in the large-jc 
gluons, while the suppression toward small pj is in line with the gluon shadowing. No addi- 
tional effects are needed to reproduce the observed spectra. It is also reassuring that this shape is 
practically independent of the fragmentation functions used in the calculation — modern sets like 
] all give equal results. 



0, 21 



Special attention should be paid to the scale-breaking effects in the data and to their good 
description by the DGLAP evolution. These effects are clearly visible e.g. in the E886 Drell-Yan 
data in Fig. ^ where the trend of diminishing nuclear effects toward larger invariant mass is 
observed. Also, from the DIS data as a function of Q^, shown in Fig. ^, the general features can be 
filtered out: At small x the 2^-slopes tend to be positive, while toward larger x the slopes gradually 
die out and become even slightly negative. 
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Figure 5: The calculated NLO scale evolution (solid black lines and error bands) of the ratio Fj^^/Fj^ and 
F2 /F2, compared with the NMC data for several fixed values of x. 

To close this section we show, in Fig. |^ a comparison of the nuclear modifications for Lead 
from all available global NLO analyses. The comparison is shown again at 2q = 1.69GeV^ and 
at = lOOGeV^. Most significant differences — that is, curves being outside our error bands — 
are found from the sea quarks and gluons. At low x, the differences shrink when the scale is 
increased, but at the high-x region notable discrepancies persist at all scales. Most of the differences 
are presumably explainable by the assumed behaviours of the fit functions, but also the choices for 
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the considered data sets (e.g. HKN and nDS do not implement the pion data), and differences in 
the definition of carry some importance. 




Figure 6: Comparison of the nuclear modifications for Lead ntQ^ — l .69GeV^ and at = lOOGeV^ from 
the NLO global DGLAP analyses HKN07 [|], nDS lU and this work, EPS09. 



4. Hadrons in the forward direction in d+Au at RHIC? 

In our previous article [^, we studied the case of a very strong gluon shadowing, motivated 
by the suppression in the nuclear modification /?dAu for the negatively-charged hadron yield in the 



forward rapidities (rj = 2.2,3.2) measured by the BRAHMS collaboration [24| in d+Au collisions 



at RHIC. We now return to this issue by applying the EPS09NLO parametrization to this specific 



process. The EPS09 predictions (with fDSS fragmentation functions []25[]) compared with the 
BRAHMS data for the absolute spectra, 

d^N^P min_bias 1 . d^N^^" n,in_bias (A^coll) ^d'o^^^ 



d^PTdy o}^^^'''^'^ d^prdy ' d^prdy a^^'^'*"'^ d^pjdy ' 

are shown in Fig. ^. In the rj =2.2 bin, the measured p+p and d+Au spectra are both in good 
agreement with the NLO pQCD. However, in the most forward rj = 3.2 bin in the p+p case, there 
is a systematic and significant discrepancy between the measured and computed px spectra. This 
observation casts doubts on any conclusion made from the nuclear modification RdAu alone — 
clearly, one should first account for the absolute baseline spectrum in p+p collisions. This is why 
we have not included this BRAHMS data set in our global analysis, either. 

5. Summary 

We have here outlined our NLO analysis of nPDFs. The very good agreement with the exper- 
imental data x'^/^ ~ 0.79 — especially the correct description of the scaling-violation effects — 
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Figure 7: Inclusive yield in p+p and d+Au collisions. The experimental pr-binned data [g4| are shown 
by open squares with statistical and systematical errors added in quadrature. The blue band indicates the 
90% confidence range derived from the CTEQ6.1M and EPS09 uncertainties. The calculated cross-sections 
have been averaged over the pr-bin width. 



lends support to the validity of collinear factorization in nuclear collisions. In addition to the best 
fit, we have distilled the experimental errors into the 30 nPDF error-sets, which encode the relevant 
parameter-space neighbourhood of the minimum. All these sets will be available as a computer 
routine at [13] for general use. 

Although not discussed here, we have also performed the leading-order (LO) counterpart of the 
NLO analysis as we want to provide the uncertainty tools also for the widely-used LO framework. 
The best-fit quality is very similar both in LO and NLO, but the uncertainty bands become some- 
what smaller when going to higher order. In the LO case, our error analysis indeed accommodates 
also the strong gluon shadowing suggested in our previous work EPS 08 [§]. 

In the near future, more RHIC data will become available and published, and factorization 
will be tested further. Also p+A runs at the LHC would be welcome for this purpose. Even better 
possibilities for the nPDF studies would, however, be provided by the lepton-ion coUiders like the 
planned eRHIC or LHeC. 

During the recent years, the free-proton PDF analyses have gradually shifted to a better or- 
ganized prescription for treating the heavy quarks than the zero-mass scheme employed in this 
analysis and in CTEQ6.1M. Such general-mass scheme should be especially important e.g. in 
the case of charged-current neutrino interactions. Interestingly, the CTEQ collaboration has lately 
looked also at this type of data among other neutrino-Iron measurements p^, and noticed that the 
best fit tends to point to somewhat different nuclear modifications than what have been obtained in 
the global nPDF analyses. Both the extension of the nPDF analysis to the general-mass scheme and 
a systematic investigation of the possible discrepancies between neutrino and other data, remains 
as future work. 
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